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(54) Title: SPACE-WEIGHTED COMMUNICATION PATH ESTIMATION 

(54) Titre: SONDAGE A PONDERATION SPATTALE D'UNE VOIE DE TRANSMISSION 



(57) Abstract 

The invention concerns a method for estimating a 
communication path formed by several channels. Said 
method requiring an estimation of the impulse response 

function (Ci, Ci C n ) of said channels comprises 

the following steps: acquiring a spatial statistic of said 
communication path; establishing a corrected impulse 
response function (C'i, C'2, C'n) at least by weighting 
the pulse response functions (Ci, C2, C n ) using the 
spatial statistic and an additive noise estimation (N01, N02, 
.... Non) of these channels. 

(57) Abrdgd 

L'invention conceme un proeddd de sondage d'une 
voie de transmission formde de plusieurs canaux. Ce 
proeddd qui ndcessite une estimation de la rdponse im- 

pulsionnelle (Ci, C2 C„) de ces canaux comprend les 

dtapes suivantes: acquisition d'une statistique spatiale de 
cette voie de transmission, dtablissement d'une rdponse 

impulsionnelle corrigde (C'i, C'2 C'n) au moins par 

ponddration des estimations de rdponses impulsionnelles 
(Ci, C2, .... C n ) au moyen de la statistique spatiale et 
d'une estimation du bruit additif (N01, N02, No n ) de 
ces canaux. 
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Sondage a ponderation spatiale d'une voie de transmission 
La presente invention concerne un procede de sondage 
d'une voie de transmission formee de plusieurs canaux. II 
s ' agit done d * une technique connue sous le nom de diversite 

5 de reception, technique selon laquelle un recepteur comporte 
plusieurs antennes chacune associee a un canal de 
transmission distinct. En d'autres terraes, 1 1 invention 
propose une raethode d' estimation de la reponse 
impulsionnelle de ces canaux de transmission. 

10 En effet, dans un systeme de transmission, notamment 

par ondes radio, le recepteur recoit pour chaque canal de 
transmission un signal emis par un emetteur. Le signal emis 
subit des fluctuations d' amplitude et de phase dans le canal 
de transmission, si bien que le signal regu par le recepteur 

15 ne lui est pas identique. Les fluctuations du signal sont 
essentiellement dues a ce que l'homme de metier appelle 
1 ' interference intersymbole . Cette interference peut 
provenir de la loi de modulation employee pour la 
transmission et elle est egalement due a la propagation 

2 0 multi-trajets dans le canal. 

II se trouve que le signal regu est generalement issu 
d'un grand nombre de reflexions dans le canal, les 
diff erents trajets empruntes par le signal emis conduisant 
ainsi a des retards varies au niveau du recepteur. La 

25 reponse impulsionnelle du canal represente 1' ensemble de ces 
fluctuations, auxquelles est soumis le signal emis. II 
s'agit done la de la caracteristique fondamentale 
representant les transmissions entre 1' emetteur et le 
recepteur ., 

30 La reponse impulsionnelle du canal est utilisee 

notamment par un egaliseur qui a precisement pour fonction 
de corriger 1 1 interference intersymbole dans le recepteur . 
Une methode classique pour realiser une estimation de cette 
reponse impulsionnelle consiste a disposer dans le signal 

35 emis une sequence d 1 apprentissage formee de symboles connus . 
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Cette sequence est choisie en fonction de la loi de 
modulation et de la dispersion du canal, dispersion devant 
s' entendre ici conraie le retard d'un symboie emis empruntant 
le trajet le plus long du canal par rapport a ce meme 
5 symboie empruntant le trajet le plus court. La dispersion 
est couramment exprimee comme un multiple de la duree qui 
separe deux symboles emis successifs, soit un nombre de 
"duree symboie" . 

A titre d'exemple, on citera deux techniques connues 
10 d' estimation de la reponse impulsionnelle d'un canal de 
transmission. 

La premiere technique fait appel a des sequences 
d 'apprentissage particulieres dites sequences CAZAC, pour 
1' expression anglo-saxonne "Constant Amplitude Zero 

15 Autocorrelation" . De telles sequences sont decrites dans 
1' article de A. MILEWSKI : "Periodic sequences with optimal 
properties for channel estimation and fast start-up 
equalization", IBM Journal of Research and Development, 
vol.27, N°5, Sept. 83, pages 426-431. 

20 Le systeme de radiocommunication cellulaire numerique 

GSM fait appel a des sequences d ' apprentissage TS formees de 
26 symboles notes a D a ^25 prenant la valeur +i ou -l . ces 
sequences possedent les proprietes suivantes : 

20 

25 2^ a i 2 = 16 

i = 5 

20 

E a i a i + k = 0 si 0 < |k| ^ 5 

i = 5 

En notant d la dispersion du canal qui vaut 4 dans le 
cas du GSM, 1' estimation de la reponse impulsionnelle prerid 
30 la forme d'un vecteur X a 5 composants notes x Q a X4. 

La sequence de symboles regus S correspondant a la 
sequence d ' apprentissage TS est formee elle aussi de 2 6 
symboles notes s D a S25. On suppose naturellement ici que 
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10 
15 
20 
25 



1 • emetteur et le recepteur sont parf aitement synchronises et 
dans ce cas 1 ' estimation de la reponse impulsionnelle X est 
donnee par 1 1 expression suivante : 



La technique CAZAC presente 1 1 avantage d • une grande 
simplicity de mise en oeuvre. Cependant, on remarque que 
chaque composante de la reponse impulsionnelle est etablie a 
partir de seulement 16 symboles recus. Etant donne que la 
sequence d 1 apprentissage comprend 26 symboles et que la 
dispersion du canal vaut 4, il y a de 1 ' information dans le 
signal regu qui n'est pas prise en compte et cela conduit a 
une reduction des performances par rapport a 1 • ideal 
theorique . 

La deuxieme technique connue fait appel au critere des 
moindres carres. Elle est decrite notamment dans les 
demandes de brevet FR 2 696 604 et EP 0 564 849. En matiere 
de rappel, cette technique fait appel a une matrice de 
mesure A construite a partir de la sequence d' apprentissage 
TS de longueur n. Cette matrice comprend (n-d) lignes et 
(d+1) colonnes, d representant tou jours la dispersion du 
canal. L ' element figurant a la ieme ligne et a la jieme 
colonne est le (d+i-j)ieme symbole de la sequence 
d 'apprentissage : 



20 




pour 0 <. k <, 4 



r a4 a3 &2 



a 5 a 4 a 3 a 2 a l 



a 6 a 5 a 4 a 3 a 2 



A = 



a 25 



a 21 
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La sequence d ' apprentissage est choisie telle que la 
matrice A*-a soit inversible ou l'operateur .t represente la 
transposition. C'est naturellement le cas pour les sequences 
CAZAC mais c'est egalement le cas pour d'autres sequences. 
5 Dans la sequence de symboles recus, on ne prend pas en 

compte les quatre premiers sq a s 3 car ceux-ci dependent 
egalement de symbole inconnus emis avant la sequence 
d' apprentissage, etant donne que la dispersion du canal vaut 
4 . Par un abus de langage on def inira done dorenavant le 

10 signal recu comme un vecteur S ayant pour composantes les 
symboles recus, S4, S5, S5, ... , s 2 5- 

Des lors, l 1 estimation de la reponse impulsionnelle 
prend la forme suivante : 

X = (At A) -1 A t . S 

15 Cette technique des moindres carres est un peu plus 

complexe que la precedente mais il faut noter que la matrice 
(A fc A) -1 A fc est calculee une seule fois. On remarque ici que 
chacune des composantes de 1' estimation de la reponse 
impulsionnelle X est etablie a partir de 22 symboles recus 

20 et non pas de 16 comme dans le cas de la technique CAZAC. On 
doit done s'attendre a une amelioration des performances. 

Cependant, quelle que soit la technique utilis6e, les 
reponses impulsionnelles de chacun des canaux de la voie de 
transmission sont considerees comme independantes . 

25 La presente invention a ainsi pour premier objet un 

procede de sondage d'une voie de transmission qui prend en 
compte le fait que les differentes antennes sont 
spatialement liees. 

Selon 1' invention, le procede de sondage d'une voie de 

30 transmission s' applique lorsque cette voie est form6e de 
plusieurs canaux et il necessite une estimation de la 
reponse impulsionnelle C-^ , C 2 , C n de ces canaux. Le 

procede comprend les etapes suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
35 transmission, 
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- etablissement d'une reponse impulsionnelle corrigee au 
moins par ponderation de ces estimations de reponses 
impulsionnelles au moyen de cette statistique spatiale et 
d'une estimation du bruit additif de ces canaux. 
5 Avantageusement, cette statistique spatiale correspond 

a une estimation de la correlation des canaux de 
transmission pris deux a deux. 

Suivant un mode de realisation pref erentiel , cette 
estimation de la correlation des canaux de transmission 

10 prend la forme d'une matrice de correlation spatiale ainsi 
definie que 1' element figurant a la ieme ligne et a la jieme 
colonne soit obtenu par lissage du produit de la transposee 
hermitienne de la reponse impulsionnelle estimee du ieme 
canal et de la reponse impulsionnelle estimee du jieme 

15 canal. 

selon une caracteristique additionnelle, un signal 
regu S par un premier canal correspondant a une sequence 
d 1 apprentissage emise, en notant A la matrice de mesure 
associee a cette sequence d ' apprentissage, 1' estimation du 
20 bruit additif Nqi de ce canal est obtenue par normalisation 
de l'energie du vecteur (S - AC^). 

on peut prevoir que cette normalisation soit suivie 
d 1 une etape de moyennage . 

De plus, une matrice de bruit N etant formee a partir 
25 de 1' estimation du bruit additif Nqi/ Nq2/ •••/ N 0ll des 
canaux, une matrice de ponderation spatiale G ' etant ainsi 
definie a partir de la matrice de correlation spatiale G et 
de la matrice de bruit : 

G' = G (G + N) _1 , 
30 les reponses impulsionnelles corrigees C'^, C' 2 , ... , C' n 
sont obtenues par 1 ' expression suivante : 
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Le procede de sondage d'une voie de transmission est 
done mis en oeuvre a partir de 1' estimation de la reponse 
impulsionnelle des differents canaux consideres comme autant 
5 de canaux independants . Or les erreurs d ' estimation sont 
inevitables. De maniere generale, la determination de la 
reponse impulsionnelle d'un canal de transmission unique est 
un probleme qui ne peut etre resolu de facon exacte en 
presence de bruit additif . De plus les techniques 

10 anterieures font l'hypothese implicite que cette reponse 
impulsionnelle peut prendre une forme quelconque. 

La presente invention a ainsi pour deuxieme objet une 
methode d* estimation de la reponse impulsionnelle d'un canal 
de transmission qui presente une meilleure resistance au 

15 bruit additif ou, autrement dit, qui mene a une erreur 
reduite comparee a 1 ' erreur d 1 estimation des techniques 
connues. Cette methode sera avantageusement appliquee pour 
l'un au moins des canaux formant la voie de transmission 
avant 1 ' etablissement des reponses impulsionnelles corrigees 

20 par ponderation au moyen de la statistique spatiale et d'une 
estimation du bruit additif de ces canaux. 

Selon 1' invention, cette methode necessite un signal 
recu par un canal , ce signal regu correspondant a une 
sequence d 1 apprentissage emise. La methode comprend les 

25 etapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission , 

- etablissement de 1' estimation de la reponse impulsionnelle 
de ce canal , cette estimation etant ponderee par la 

30 statistique temporelle du canal au moyen du signal regu. 
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La statistique temporelle du canal represente une 
valeur de la reponse impulsionnelle anterieure a 



introduit le fait que la reponse impulsionnelle afferente au 
5 signal recu a une valeur probablement plus proche de cette 
valeur anterieure qu'une valeur qui en serait tres eloignee. 
Ainsi, statistiquement , l'erreur d' estimation est diminuee. 

Avantageusement , cette statistique correspond a une 
estimation de la covariance de cette reponse impulsionnelle. 
10 Selon une premiere variante du procSde, celui-ci 

comprend les etapes suivantes : 

- lissage de la reponse impulsionnelle et orthonormalisation 
au moyen d'une matrice de transformation W pour obtenir 
1 1 estimation de la covariance qui prend alors la forme d • une 

15 matrice L 1 , 

- recherche des vecteurs propres v^ 1 et valeurs propres Xj_ 1 
associees de cette matrice L 1 , 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir du signal regu et application de cette 

20 matrice de transformation W pour former un vecteur X' , 

1 ' etablissement de 1' estimation ponderee Xp etant ainsi 
realised s 



ou Nq est un nombre r6el strictement positif representant le 
25 bruit additif du canal. 

On peut ici prevoir que le bruit additif soit choisi 
egal a la plus petite des valeurs propres X.^ • . 

Lorsqu'un sous-ensemble de ces valeurs propres X^' 
presente une contribution inferieure a un seuil 
30 predetermine, on peut egalement prevoir que chacune de ces 
valeurs propres soit f orcee a la valeur du bruit additif . 



Selon une seconde variante du procede, 1' estimation de 
la covariance prenant la forme d'une matrice R, en notant A 



1 'acquisition du signal regu. La susdite pond^ration 




La complexite s 1 en trouve reduite d'autant. 
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la matrice de mesure associee a la sequence d ' apprentissage , 
1 ' etablissement de 1' estimation ponderee est ainsi realise : 

Xp = (A fc A + NqR -1 ) -1 At.S 
ou N 0 est un reel strictement positif representant le bruit 
5 additif du canal . 

Par ailleurs, il est possible d'effectuer une etape 
d ' orthonormalisation de la matrice R au moyen d'une matrice 
de transformation W pour obtenir une nouvelle matrice R • , 
1' estimation ponderee prenant alors la nouvelle forme 
10 suivante : 

Xp = W fc (I + NqR"- 1 )- 1 W fc A't.S 
ou la matrice A* est egale au produit de la matrice de 
transformation W et de la matrice de mesure A. 

Avantageusement, 1' expression (I + NqR' -1 ) -1 est 
15 calculee au moyen du lemme d 1 inversion matricielle. 

La presence invention apparaitra maintenant de maniere 
plus detaillee dans le cadre de la description qui suit oil 
sont proposes des exemples de mise en oeuvre a titre 
illustratif, ceci en reference aux figures annexees qui 
2 0 representent ; 

la figure 1 , un diagramme identif iant les 
principales etapes d'un mode de mise en oeuvre du procede 
de sondage d'une voie de transmission selon 1 ' invention , 

- la figure 2, une premiere variante d' estimation de 
25 la reponse impulsionnelle d'un canal selon 1' invention, 

- la figure 3, une deuxieme variante d 1 estimation de 
la reponse impulsionnelle d'un canal selon l 1 invention. 

L 1 invention sera presentee dans son application au GSM 
car ce systeme a le merite d 1 etre bien connu de 1 ' homme du 
30 metier. II s'agit done la d'une presentation adoptee dans un 
souci de clarte mais il ne faut y voir en aucun cas une 
limitation de 1' invention a ce seul systeme. 

En reference a la figure 1, le procede de sondage 
d'une voie de transmission s' applique lorsque cette voie 
35 comporte au moins deux canaux de transmission, n canaux dans 
le cas general. Les canaux correspondent chacun a une 
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antenne distincte. On considere done que chaque canal a fait 
l'objet d'une estimation de sa reponse impulsionnelle 
respective C lr C 2 , . .., C n selon l'une quelconque des 
techniques disponibles. 
5 En premier lieu le procede prevoit 1 ' acquisition d'une 

statistique spatiale de cette voie de transmission. Par 
statistique spatiale on entend un ensemble de donnees 
refletant le comportement de cette voie sur une periode 
predetermines que l'on conviendra d'appeler periode de 

10 correlation. En effet, les differentes antennes etant fixes, 
les signaux recus sur ces antennes presentent un certain 
degre de correlation. L' invention entend precisement tirer 
parti de ce fait pour ameliorer la qualite d' estimation de 
la reponse impulsionnelle de l'un des canaux au moins. A 

15 titre d'exemple, une telle statistique peut etre obtenue au 
moyen d'une matrice de correlation spatiale G : 



''vecj^) V(C^C 2 ) .......... V(C^C n ) >l 

VCC^) V(C^C 2 ) v ( C 2 C n ) 



G = 



4 • O • » • 



v(cjjc 1 ) v(c^c- 2 ) 



v(chc n ) 



ou l'operateur . h represente la transposition hermitienne. 
20 Cette matrice G carr^e de dimension (n,n) peut ainsi 

etre representee de maniere generique par 1 • element g^ j 
figurant a la ieme ligne et a la jieme colonne ; 
g. . = V(C^C. ) 

25 au moyen des reponses impulsionnelles estimees Ci , Ch des 



Cet element gij est obtenu par lissage du produit C^C^ 

Ci , C- 



au moyen des reponses impulsionnelles estimees <^ Xr 
ieme et jieme canaux obtenues durant la periode de 
correlation. Ce lissage est une estimation de la correlation 
de ces deux canaux. 
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On entend ici lissage dans un sens tres general , 
c 1 est-a-dire toute operation permettant de lisser ou de 
moyenner le produit c|?C. sur la periode de correlation. 

Un premier exemple de lissage consiste a faire la 
5 moyenne de ce produit sur cette periode supposee comprendre 
e estimations successives de chacune des reponses 
impulsionnelles C^, C2, . . . r C n : 

e 

v(c^c. ) = -Yck 

1 

10 Un deuxieme exemple de lissage consiste a actualiser, 

a la pieme estimation recue pour chacun des ieme et jieme 
canaux, 1' expression de lissage v p _ i( c i c j ) obtenue a la 

( p - 1 )ieme estimation au moyen d * un coefficient 
multiplicatif a, ce coefficient compris entre 0 et l etant 
15 souvent connu sous le nom de facteur d'oubli de lissage 1 

V p (C i C j } = aC i C j + ^- a > v p - l (C i C j ) 
L 1 initialisation peut se faire par tous moyens, 
notamment au moyen de la premiere estimation obtenue ou bien 
par une moyenne obtenue comme dans le premier exemple pour 
2 0 les premieres estimations recues. 

Le procede de sondage propose ensuite de realiser une 
estimation du bruit additif Nqi, Nq2/ Nqh present dans 

chacun des canaux au moyen des estimations des reponses 
impulsionnelles respectives Ci, C2, C n de ces canaux. 

2 5 II sera propose differentes solutions pour mener a 

bien cette estimation du bruit et on traitera le cas d'un 
seul canal, le premier par exemple, etant donne que les 
memes solutions s'appliquent a chacun des canaux. 

Une solution simple consiste a affecter Nqi d'une 
30 valeur predetermined qui reflete un seuil en dessous duquel 
il est peu probable que le bruit additif puisse descendre. 
Cette valeur pourrait etre determinee par une mesure de 
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rapport signal a bruit, ou par les performances du 
recepteur, ceci a titre d'exemple. 

Par ailleurs, 1 1 estimation du bruit Nqi peut etre 
obtenue au moyen de 1 ' estimation de la reponse 
5 impulsionnelle du premier canal et du signal recu 

correspondant S. Cette estimation du bruit a done pu etre 
realisee avant la mise en oeuvre du procede de sondage a 
ponderation spatiale selon quelque technique que ce soit . 
Toutefois, si ce n'etait pas le cas, on propose ici une 
10 methode appropriee lorsque 1' estimation de la reponse 
impulsionnelle est acquise selon la technique des 

moindres carres . 

II s'agit done d'estimer directement le bruit additif 
a partir du signal recu S et de la matrice de mesure A. En 
15 effet, en notant U 1 le vecteur bruit affectant le signal 
recu, il vient que : 

S = AC! + N X 

Compte tenu du fait que les vecteur s S et N-l ont 2 2 
composantes, 1' estimation du bruit Nqi prend la forme 
20 suivante ; 

N 0 1 = ("4> (S " ACi)h ( S - AC!) 

Naturellement cette estimation du bruit additif Nqi 
peut etre moyennee ou lissee. 

A partir des estimations du bruit additif des 
25 differents canaux N 0 i, Nq2^ •••/ N 0n' on construit une 
matrice de bruit N dont 1 ' element figurant a la ieme ligne 
et a la jieme colonne est note njj . 

Lorsque i est egal a j , cet element n^j prend la 
valeur du bruit estime Nqi- 
30 Lorsque i et j sont differents, 1' element n^j est nul . 

Par ailleurs, toujours dans un souci de 
simplification, il est possible de calculer la valeur 
moyenne N 0 des bruits estimes Nqi, Nq2/ • • • > N 0n et de 
forcer chacun des elements n^j de la diagonale de cette 
35 matrice N a cette valeur moyenne Nq. En notant I la matrice 
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identity, il s'ensuit que la matrice de bruit N prend la 
forme suivante : 

N = N 0 I 

La matrice de correlation G et la matrice de bruit N 
5 permettent de definir une nouvelle matrice, la matrice de 
ponderation spatiale G 1 : 

G ' = G(G + N) _1 
Les reponses impulsionnelles corrigees C'^, C'2/ 
C* n au moyen de cette matrice de ponderation spatiale G' 
10 sont ainsi definies : 
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ou. l'operateur represente tou jours la transposition. 

L'une au moins des reponses impulsionnelles corrigees 
15 C*i est utilisee a la place de 1' estimation de la reponse 

impulsionnelle dans le recepteur. 

A titre optionnel, 1' invention permet egalement 

d'ameliorer 1' estimation de la reponse impulsionnelle C]_, 

C2, . = C n de chacun des canaux, ceci avant la mise en 
20 oeuvre du procede de sondage a ponderation spatiale d'une 

voie de transmission, c'est-a-dire avant 1 1 etablissement de 

la ou des reponses impulsionnelles corrigees C'^, C'2/ 

C ' n . 

On prevoit ainsi 1 'acquisition d'une statistique 
25 temporelle pour l'un au moins de ces canaux, le premier par 
exemple. Par statistique temporelle, on entend un ensemble 
de donnees refletant le comportement de ce canal considere 
independamment des autres canaux sur une periode 
predeterminee que 1 ' on conviendra d ' appeler periode 
3 0 d' analyse. II s'agit done d'une representation du 
comportement moyen du canal pendant la periode d' analyse. 
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Cette statistique peut etre etablie par quelque moyen que ce 
soit et en quelque lieu que ce soit. En effet, 
1 ' etablissement de la statistique peut prendre place dans un 
equipement quelconque du reseau de radiocommunication . Ce 
5 qui importe c'est que le recepteur puisse acquerir cette 
statistique . 

A titre d'exemple, une telle statistique temporelle 
peut etre obtenue de la maniere suivante. 

Au moyen du signal recu S pendant la periode 
10 d' analyse, on calcule une estimation X de la reponse 
impulsionnelle selon une methode connue. 

Si l'on retient la technique des moindres carres cette 
estimation X vaut : 

X = (A fc A) -1 At.S (1) 
15 On rappelle ici que l'emetteur et le recepteur sont 

supposes synchronises a mieux que un demi-symbole pres , 
auquel cas le signal de reception est forme par le vecteur S 
dont les composantes sont les symboles regus S4 a S25 
synchrones des symboles a4 a a25 de la sequence 
2 0 d 1 apprentissage TS . Si une telle synchronisation n * etait pas 
acquise, plusieurs solutions sont disponibles pour 
1 1 acquerir et on en citera deux a titre d 1 exemple . 

La premiere solution consiste a decaler le signal regu 
en avance ou en retard de j periodes symboles de sorte que 
25 Sjt = (s 4 _j, s 5 _ j , s 6 _j / ... s 2 5-j ) . 

On calcule alors 1' estimation Xj pour chacun des 
vecteurs Sj et on retient la valeur pour laquelle Xj n °Xj 
est maximal . Cette valeur donne la synchronisation 

escomptee et il suffit de remplacer le vecteur S dans 
30 l 1 equation (1) par le vecteur Sj^. 

La deuxieme solution consiste a augmenter 
artif iciellement la dispersion d du canal d'une quantite "2q 
predeterminee . On peut alors definir une matrice de mesure 
modifiee A m comprenant (n-d-2q) lignes et (d+2q+l) colonnes. 
35 En reprenant les valeurs 26 et 4 respectivement pour n 
et d : 
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Am = 



( 



a 4 + 2q 
a 5 + 2q 
a 6 + 2q 
a 7 + 2q 



o e • o 



a 4 a 3 a2 ai ao 
a 5 a 4 a 3 a 2 a l 
a6 as a4 a3 a2 



00004000 



a 25 



a 21 - 2q 



II faut 

alors reduire le nombre de composantes du signal recu S de 
cette meme quantite 2q et on peut convenir de retenir le 
5 vecteur S m modif ie ; 

s m^ = ( s *4+q/ s ' 5+q/ • • • s '25-q) 
On obtient ainsi line estimation modifiee X m % 

% = (A m fc Am) Amt-Sm 

Cette estimation modifiee X m comporte d+2q+l 

10 composantes : 

Xjn^ = (x_q, o . . Xq / / . . . , X4/ = » . x^^g) 

On recherche alors la valeur de j comprise entre -q 

et +q qui maximise 1' expression suivante : 
4 

2 xt + k xj +k 
k = o 

15 ou l'operateur .* represente la conjugaison complexe. 

Cette valeur determine 1' estimation X de la reponse 
impulsionnelle pour une dispersion d=4 : 

xt = ( x jM' x jM+l/ ••• x jM+4) 
La synchronisation s ' en deduit immediatement en 
20 appliquant le decalage au signal recu S. 

On peut des lors appliquer a nouveau 1 1 equation ( 1 ) . 
On construit maintenant une matrice de lissage L par 
lissage des differentes estimations X obtenues pendant la 
periode d' analyse pour obtenir une estimation de la 
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covariance associee a cette reponse impulsionnelle. On 

entend ici lissage dans un sens tres general, C 1 est-a-dire 

toute operation permettant de lisser ou de moyenner la 

reponse impulsionnelle sur la periode d' analyse. On obtient 

5 ainsi une representation statistique du comportement du 

canal de transmission. Ce lissage peut etre obtenu notamment 

par l'une des deux methodes proposees ci-dessus, 

1 ' expression la plus simple de cette matrice de lissage L 

etant la suivante : 

i m 

10 L = -V XX h 

ra *Y 

ou m correspond au nombre de sequences d 1 apprentissage sur 
lequel est calcule le lissage. 

On admet ici que cette matrice de lissage peut etre 
approchee par 1 • equation suivante : 
15 L ~ (AtA) -1 N 0 + R (2) 

ou Nq represente tou jours le bruit present dans le 
canal de transmission ou bruit additif et ou R est une 
matrice que l'on a coutume d'appeler statistique a priori du 
canal car elle represente le comportement du canal 
20 abstraction faite du bruit. 

On admet egalement que la matrice de mesure A est bien 
conditionnee, c' est-a-dire que les valeurs propres de la 
matrice A^A. sont tres proches les unes des autres . Dans ce 
cas, il est interessant de proceder a 1 •orthonormalisation 
25 des vecteurs constitues par les lignes de la matrice de 
mesure A, mais il ne faut pas voir la une limitation de 
1 1 invention . 

Pour ce faire, on emploie une matrice de 
transformation W telle que : 
30 A' = AW et A ,fc A' = I 

ou I represente la matrice identite. 

En notant L' la matrice ainsi definie : 
L = WL'W fc , 

il vient que 1' equation (2) peut maintenant s'ecrire : 
35 L' ~ N 0 I + R 1 (3) 
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Selon une premiere variante representee a la figure 2 , 
on remarque que les vecteurs propres 1 de L' et de R' 
sont identiques tandis que les valeurs propres X^' de L' et 
A,i de R 1 sont decalees de Nq. Soit en prenant toujours la 
5 meme valeur de 4 pour la dispersion du canal, pour tout i 
cornpris entre 0 et 4 : 
Vi« = V ± 
Xi« = Xi + N 0 

II apparait ainsi que la determination des vecteurs 
10 propres et valeurs propres de R' et celle de L' sont 
identiques sous reserve que Nq soit connu. 

L'etape d 1 estimation du bruit sera decrite plus loin 
pour rendre 1 ' expose plus clair, bien que cette 6tape 
precede celle qui va maintenant etre explicitee. 
15 II convient done maintenant de rechercher les couples 

valeur-propre/vecteur-propre pour l'une ou 1' autre des 
matrices L' ou R 1 . Cette etape ne sera pas plus detaillee 
car bien connue de l'homme du metier. Par ailleurs, il va 
sans dire que l'on peut annuler les valeurs propres dont la 
20 contribution est jugee non significative. Par exemple, si 
ces valeurs propres sont classees par ordre decroissant, on 
supprime les dernieres qui sont telles que leur somme soit 
inferieure a un seuil predetermine. 

L'etape suivante consiste a estimer la reponse 
25 impulsionnelle instantanee X etablie selon l'une quelconque 
des techniques connues a partir du signal regu correspondant 
a la derniere sequence d' apprentissage regue. En notant 
X = WX 1 , cette derniere estimation est ponderee par la 
methode suivante pour obtenir une ponderation temporelle Xp 
3 0 de la reponse impulsionnelle instantanee : 
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Pour obtenir la ponderation Xp il faut done estimer le 
bruit additif Nq. 

Cette estimation du bruit peut etre effectuee 
notamment selon 1 ' une des methodes exposees plus haut . 
5 Une autre possibility consiste a considerer que la 

derniere valeur propre, (la plus faible) de la matrice de 
lissage L est egale a Nq . 

A-4 1 = Nq ou X4 = 0 . 
Quelle que soit la methode retenue, la ponderation 
10 temporelle Xp de 1 1 estimation de la reponse impulsionnelle 
instantanee peut alors etre realisee comme mentionne ci- 
dessus . 

Selon une seconde variante representee a la figure 3, 

on etablit 1' estimation ponderee Xp directement comme suit : 

15 Xp = (A fc A + NoR -1 ) - ^^ 

ou bien en reprenant la matrice de transformation W 

definie ci-dessus t 

Xp = W (I + NqR' -1 ) -1 W fc A !t .S (4) 

Conformement a l 1 equation (3) : 

20 R' = L 1 - N 0 I 

II convient la aussi d* estimer le bruit additif Nq. 

Une solution avantageuse pour obtenir la ponderation 

temporelle Xp consiste a adopter la methode suivante. 

On divise la matrice R ' par Nq : 
r * 

25 B = — 

N 0 

II s'ensuit que ; 

I + NqR'- 1 = I + B _1 
On utilise le lemme d' inversion matricielle pour 
calculer la matrice de pond6ration P = (I + B -1 )" 1 . 
30 Ainsi, en notant e^ les vecteurs canoniques, on 

procede a 1 1 iteration suivante : 

- initialisation : 

P = B 

- pour i variant de 0 a d (4 dans le cas present) ; 
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Pe (pe ) 

P = P r— = 

1 + e?Pe. 
x 1 

P etant connu, il reste a etablir la ponderation Xp 
selon 1 1 equation ( 4 ) . 

On remarquera que la matrice de ponderation P n'est 
5 pas necessairement calculee a l 1 occasion de chaque nouvelle 
sequence d ' apprentissage emise. Elle peut etre calculee a un 
rythme plus lent car elle varie sensiblement au meme rythme 
que R" et done plus lenteraent que le signal regu S. 

On remarquera egalement que 1 1 estimation ponderee est 
10 realisee sans avoir recours a la reponse impulsionnelle 
instantanee. Elle est produite directement a partir du 
signal recu S . 

Quelle que soit la variante retenue, e'est done 
1 1 estimation ponderee Xp qui sera avantageusement utilisee 
15 comme estimation de la reponse impulsionnelle pour la 

mise en oeuvre du procede de sondage a ponderation spatiale 
d'une voie de transmission, e'est-a-dire pour produire une 
ou plusieurs reponses impulsionnelles corrigees. 
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REVEMDIC&TIOMS 

1) Procede de sondage d'une voie de transmission 
formee de plusieurs canaux, ce procede necessitant une 
estimation de la reponse impulsionnelle (C^, C21 • .. /C n ) de 

5 ces canaux, caracterise en ce qu'il comprend les etapes 
suivantes 5 

- acquisition d'une statistique spatiale de cette voie de 
transmission , 

- etablissement d'une reponse impulsionnelle corrigee (c'i, 
10 C r 2/ f c 'n) au moins P ar ponderation desdites estimations 

de reponses impulsionnelles (C^, C2/ ... ,C n ) au moyen de 
ladite statistique spatiale et d'une estimation du bruit 
additif (N 0 i, Nq2/ N 0n) desdits canaux. 

2) Procede selon la revendication 1, caracterise en ce 
15 que ladite statistique spatiale correspond a une estimation 

de la correlation desdits canaux de transmission pris deux a 
deux. 

3) Procede selon la revendication 2, caracterise en 
ce que cette estimation de la correlation des canaux de 

20 transmission prend la forme d'une matrice de correlation 
spatiale (G) ainsi definie que 1' element (g±j) figurant a la 
ieme ligne et a la jieme colonne soit obtenu par lissage du 
produit ( C^C.. ) de la transposee hermitienne de la reponse 

impulsionnelle estimee (C^) du ieme canal et de la reponse 
25 impulsionnelle estimee (Cj) du jieme canal. 

4) Procede selon la revendication 3 caracterise en 
ce que, un signal regu S par un canal correspondant a une 
sequence d'apprentissage emise, en notant A la matrice de 
mesure associee a ladite sequence d'apprentissage, 

30 1' estimation du bruit additif (Nqi) de ce canal est obtenue 
par normalisation de 1 ' energie du vecteur (S - ACi). 

5) Procede selon la revendication 4, caracterise en ce 
que ladite normalisation est suivie d'une etape de 
moyennage . 

35 6) Procede selon l'une quelconque des 

revendications 3 a 5 caracterise en ce que, une matrice de 
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bruit (N) etant formee a partir de 1* estimation du bruit 
additif (N 01 , N 02r — , N 0n ) des canaux, une matrice de 
ponderation spatiale (G 1 ) etant ainsi definie a partir de 
ladite matrice de correlation spatiale (G) et de ladite 
5 matrice de bruit : 

G" = G (G + N) -1 , 
lesdites reponses impulsionnelles corrigees (C'i, C'2, • • . / 
C' n ) sont obtenues par 1" expression suivante : 
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10 7) Procede selon l'une quelconque des 

revendications 1 a 6 caracterise en ce que, avant 
1 1 etablissement de ladite reponse impulsionnelle corrig#e 
(C 1 !) d'un canal, le signal regu (S) par ce canal 
correspondant a une sequence d ' apprentissage emise, il 

15 comprend les etapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission , 

- etablissement de 1' estimation de la reponse impulsionnelle 
(Ci) dudit canal, cette estimation (Xp) etant ponderee par 

20 la dite statistique temporelle du canal au moyen dudit 
signal recu (S). 

8 ) Procede selon la revendication 7 , caracterise en 
ce que ladite statistique temporelle correspond a une 
estimation de la covariance de ladite reponse 

25 impulsionnelle. 

9 ) Procede selon la revendication 8 caracterise en 
ce qu'il comprend les etapes suivantes : 

- lissage de ladite reponse impulsionnelle et 
orthonormalisation au moyen d'une matrice de transformation 

3 0 W pour obtenir ladite estimation de la covariance qui prend 
alors la forme d'une matrice L' , 
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- recherche des vecteurs propres (Vj_') et valeurs propres 
(A-i') associees de cette matrice L', 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir dudit signal recu (S) et application de cette 
matrice de transformation W pour former un vecteur X' , 

1 ' etablissement de ladite estimation ponderee (Xp) etant 



ainsi realisee 



Xp = 



V A i J 



W V. 11 
1 



ou No est un nombre reel strictement positif representant le 
10 bruit additif dudit canal. 

10) Procede selon la revendication 9, caracterise en 
ce que ledit bruit additif (Nq) est choisi egal a la plus 
petite desdites valeurs propres (A,^*)„ 

11) Procede selon l'une quelconque des revendications 9 
15 ou 10, caracterise en ce qu'un sous-ensemble desdites 

valeurs propres (^j.') presentant une contribution inferieure 
a un seuil predetermine, chacune de ces valeurs propres est 
forcee a la valeur dudit bruit additif (Nq). 

12) Procede selon la revendication 8 caracterise en 
20 ce que ladite estimation de la covariance prenant la forme 

d'une matrice R, en notant A la matrice de mesure associee a 
ladite sequence d'apprentissage, 1 ' 6tablissement de ladite 
estimation ponderee (Xp) est ainsi realise : 
Xp = (A fc A + NqR -1 ) -1 At.S 
25 ou Nq est un reel strictement positif representant le bruit 
additif dudit canal . 

13) Procede selon la revendication 12, caracterise 
en ce qu'il comprend une etape d • orthonormalisation de 
ladite matrice R au moyen d'une matrice de transformation W 

3 0 pour obtenir une nouvelle matrice R 1 , cette estimation 
ponderee prenant alors la nouvelle forme suivante : 

Xp = W fc (I + NqR' -1 ) -1 W fc A't.S 
ou la matrice A 1 est egale au produit de cette matrice de 
transformation W et de ladite matrice de mesure Ac 
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14) Procede selon la revendication 13, caracterise 
en ce que 1' expression (I+NqR 1-1 ) -1 est calculee au moyen du 
lemme d' inversion matricielle . 
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